Ground State Geometry of Binary Condensates in Axisymmetric Traps 



S. Gautam and D. Angom 
Physical Research Laboratory, 
Navarangpura, Ahmedabad - 380 009 

(Dated: November 24, 2009) 

We show that the ground state interface geometry of binary condensates in the phase separated 
regime undergoes a smooth transition from planar to ellipsoidal to cylindrical geometry. This 
occurs for condensates with repulsive interactions as the trapping potential is changed from prolate 
to oblate. The correct ground state geometry emerges when the interface energy is included in the 
energy minimization. Where as energy minimization based on Thomas-Fermi approximation gives 
incorrect geometry. 
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I. INTRODUCTION 

Two species Bose-Einstein condensate (TBEC), con- 
sisting of two different hyperfine spin sates of 87 Rb, was 
first observed by Myatt et al pQ. Since then, TBECs of 
different atomic species ( 41 K and 87 Rb) [5] and of dif- 
ferent isotopes of the same atomic species [5] have been 
experimentally realized. This has led to several experi- 
mental and theoretical investigations on different aspects 
of TBECs. The remarkable feature of TBECs which 
is absent in a single component Bose-Einstein conden- 
sates (BECs) is the phenomenon of phase separation. In 
Thomas- Fermi approximation (TEA), the phase separa- 
tion occurs when all the inter atomic interactions are 
repulsive and the inter species repulsion exceeds the ge- 
ometric mean of the intra species repulsive interactions. 

Depending upon the properties of the condensates 
and trapping potential parameters, the ground state of 
TBECs assumes a configuration which minimizes the to- 
tal energy. The structure of the ground state plays an 
important role in dynamical phenomena like Rayleigh- 
Taylor [4, 5 instability, Kelvin- Hclmholtz instability |6J, 
pattern formation at the interface [7] etc. It was recently 
demonstrated that quantum RaylcighTaylor instability 
can be observed in a very controlled way with TBECs in 
cigar shaped traps [4]. 

In the phase separated regime, the interface energy of 
the two component species defines the geometry of the 
ground state. In a previous work, the ground state ge- 
ometry of TBECs was examined within the TFA j8] that 
is without the interface energy. In later works, [HI EH EQ 
the contribution from the interface energy was incorpo- 
rated. From these it is observed that the analytic approx- 
imations for interface energy are not sufficient enough to 
explain the experimental results of strongly segregated 
ground states [3] . A recent work [T^] reported a more ac- 
curate determination of the interface energy. It explains 
the stationary state geometry of the strongly segregated 
TBECs more precisely. 

In this paper we provide a semi-analytic scheme to 
determine the stationary state structure of TBEC in ax- 
isymmteric traps. For this, we follow the ansatz adopted 



in Ref. (T3] i.e. to minimize the total energy of TBEC 
with fixed number of particles of each species in TFA. 
In section |n] of the manuscript, we identify three geome- 
tries which a TBEC can assume depending on the trap- 
ping potential parameters. Based on TFA, we provide a 
semi-analytic scheme to determine the stationary state 
parameters of the ground state for each of these three 
geometries. In section |III| we analyze the crucial role 
played by the interface energy in determining the ground 
state structure of the TBEC. 



II. TBEC IN AXISYMMETRIC TRAPS 

We consider TBEC in axisymmetric trapping poten- 
tials 

Vl (r,z) = ^f(ay + \i Z i), (1) 

where i = 1,2 is the species index, and ai and are 
the anisotropy parameters. In the mean field approxima- 
tion, the stationary state solution of binary condensate is 
described by a set of coupled Gross-Pitaevskii equations 

-^V 2 + Vi(r,z) + U u \A(r,z)\ 2 + U tJ \^(r,z)\ 2 
2m,i 

i>i(r,z) = [Hil>i(r,z), (2) 

here i and j = 3 — i are species indices; Uu = 4irh 2 ai/rrii 
with TO, as mass and as s-wave scattering length, is 
the intra-species interaction; Uij — 2irH 2 a,ij / rriij with 
rriij = rriimj /{mi + rrij) as reduced mass and ay as inter- 
species scattering length, is the inter-species interaction 
term and Hi is the chemical potential of the ith species. 

When the number of atoms are large, the TFA is appli- 
cable to obtain the stationary state solutions of Eq.Q. 
In this limit the kinetic energy is neglected in compar- 
ison to interaction energy. We consider the interaction 
parameter U±2 > VUrJ-h2, such that the two compo- 
nents are phase separated. That is, the two components 
occupy different regions of the trapping potentials. Ne- 
glecting the overlap between the species, stationary state 
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solutions within TFA are 



\ipi(r,z)\ 



. , _ fj,j - Vj(r, z) 



(3) 



where fii is fixed by the number of atoms of the corre- 
sponding species. 

The total energy of the TBEC in the phase separated 
regime is 



E = / dV[V 1 (r,z)\Tp 1 (r,z)\ 2 + V 2 (r,z)\i> 2 (r,z)\ 2 + 



(4) 



Depending upon the anisotropy parameters, the TBEC 
can have three distinct spatial distributions in axisym- 
metric traps. The distinguishing feature of these struc- 
tures is the geometry of the interface, which can be pla- 
nar, cylindrical or ellipsoidal. The smooth transition of 
interface geometry, for the TBEC of 85 Rb- 87 Rb mixture, 
from planar to ellipsoidal and finally to cylindrical is 
shown in Fig. [T] These features are most prominent 
when the TBEC is strongly segregated and for the de- 
tailed examination of our scheme we choose 85 Rb- 87 Rb 
experiments of Papp et al [3] . Where two of the geome- 
tries, planar and ellipsoidal, were observed. 



A. Planar interface 

It has been observed experimentally [3] that in cigar 
shaped traps (Aj <C a«) the TBEC assumes a sandwich 
structure with planar interface between the two species. 
In this structure the phase separation occurs along the 
axial direction, with the weakly interacting component 
sandwiched by the strongly interacting one. There are 
two realizations of this: coincident and shifted trapping 
potentials. 



1. Coincident trap centers 

An idealized choice of V* is with coincident centers. If 
z = ±Li are the locations of the planes separating the 
two components and L 2 , the axial size of binary conden- 
sate. Then the problem of determining the structure of 
the TBEC is equivalent to calculating L\. If Ni and Ri 
are the number of atoms and radial size of the ith species 
respectively, then 



Ni=2ir rdr dz\ip t (r, z)\ 2 . 

JO J-Li 



(5) 



FIG. 1: The false color coded images of |V>i(r, z)\, for TBEC 
consisting of 85 Rb (inner component) and 87 Rb (outer compo- 
nent), on r—z plane with vertical and the horizontal directions 
corresponding to radial and axial coordinate respectively. The 
wave functions are obtained by numerically solving Eq.|2]| us- 



ing a a 



and ai2 taken from Ref. 



referred to as set a in 



the text. The images correspond to different values of Ai but 
same Ni = 50,000. In the first row, starting from left (Ai, 
A 2 ) are (0.020, 0.022), (0.50, 0.50) and (0.85,0.85). While 
the second row corresponds to (1.0, 1.0), (1.5, 1.5) and (50.0, 
50.0). As is evident the interface geometry changes continu- 
ously from planer to ellipsoidal to cylindrical. 



From Eq.(|3|, we get 
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Similarly, the total energy in Eq.Q is 

47T 



E = 



1680w 2 m 2 U 22 a\ 

\\2^L\m\\\pi 2 ~ 112w 4 LfmlA^2 



UOuj 2 L\m 2 \\[L\ 
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(7) 



Here, L\ is determined by variational minimization of E 
with L\ as the variational parameter and constraints that 
jLti, and /i2 satisfy Eq.([6]) for fixed Ni. In the constraint 
equations, we invert the expression of Ni and obtain 
as a function of L\. However, inverting 7V2 to calculate 
fi 2 is nontrivial. Hence, we implement the minimization 
numerically. 

As mentioned earlier we consider the TBEC of 85 Rb- 
87 Rb with Ni = 50,000. The scattering lengths a x = 
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FIG. 2: Surface and contour plot of the variation of E as 
function of Zi and L\ . The minima of E occurs at l\ — 39.5a osc 
and Li = 26.0a osc . These are close to li = 40.15a osc and L\ = 
27.95a osc obtained from the numerical solution of coupled GP 
equations. 



51et , a 2 — 99a Q and a 12 — 214et are from the exper- 
imental results of Wicman and collaborators [3 . Like 
wise the anisotropy parameters and trap frequency are 
at = 1, Ai = 2.9/130, A 2 = 2.6/130 and uj = 130Hz 
respectively. From here on this choice of parameters is 
referred as the set a. For these parameters, the minima 
of E occurs at 32.3a osc . This is in very good agreement 
with the value of 33.8a osc obtained from the numerical 
solution P3] of Eq.Q. The contour plots showing the 
absolute value of wave functions of 85 Rb and 87 Rb, ob- 
tained by numerically solving Eq.Q, are shown in first 
image from left in upper panel of Fig. [l] 



2. Separated trap centers 

In the experimental realizations, the gravitational po- 
tential and tilts in the external field configurations tends 
to separate the minima of the effective potentials. Nor- 
mally, in cigar shaped traps, the tilt angle is small and 
separation is effectively along the axial direction. Then 
potentials with separation zq are 



B. Ellipsoidal Interface 

As the anisotropy parameter A is increased, beyond 
a critical value A a the interface geometry change from 
planer to ellipsoidal. Where one species envelopes the 
other. This is the preferred interface geometry, for the 
phase separated TBEC in axisymmteric traps, without 
interface energy. Consider trapping potentials with coin- 
cident centers, if Ri and Li are equatorial (along radial 
direction) and polar (along axial direction) radii of the 
ith species respectively, then 



Ni = 2n rdr dz\ipi(r, z)\ 2 

JO J-Li 

From Eq.(|3]) and Eq.Q, we get 

-27ri?iai(3w 2 mii?fa? - 10iJ?Mi) 
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V 2 {r,z) 



UI2L0 



a y + \ 2 2 { z - z y 



(8) 



-WR iai \ lt l2 + 2 3 3/2 ' 

E ^ 210^ 2 J llt / 22 a|Af (-^ 6 mlm 2 RlU 22 alalXi 
+Lu 6 m 2 i RlU 11 a 1 al(8a 2 i \i + ±a\a\\\\\ + Saf\ 2 ) 
+160V2U 11 X 5 ll i 3 2 J / 2 

+ UOuj 2 m 2 R 3 1 a 1 ajX 4 1 (U 22 fi 2 1 - I7 n ^)) . (12) 

In TFA, the profile of density \ipi(r, z)\ 2 has the same el- 
lipticity e as that of the trapping potential, which is a 
function A. Then, Li is a.iRi/Xi, further more in TFA fi 2 
constrains the value of R 2 . These reduce the variation 
parameter to only R\. The energy E is then minimized 
numerically to find the equilibrium geometry of the phase 
separated TBEC. To examine the scheme consider the 
85 Rb- 87 Rb mixture with parameter set a and coincident 
trapping potentials, however take Ai as 1.5. Then the 
equilibrium geometry is ellipsoidal with an equitorial ra- 
dius (i?i ) of 3.72a OS c- This is in very good agreement 
with value 3.75a osc obtained from the numerical solution 
of GP equations. The contour plots showing the absolute 
value of wave functions of 85 Rb and 87 Rb, obtained by 
numerically solving Eq.Q, are shown in second image 
from left in lower panel of Fig. [TJ 



Due to the lack of axial symmetry, z = —l\ and z = 
L\ are the two planes separating the two components. 
These, l\ and L\, are the parameters to minimize E. 
Like in the previous subsection N\, N 2 and E can be 
evaluated and are presented in the appendix. For the 
parameter set a and zq — 3.4/im, the minima of E occurs 
when lx and L\ are 39.5a osc and 26.0a osc respectively. 
The overall trend of E is shown in Fig. [2j 



C. Cylindrical interface 

On further increase of A, beyond another critical value 
Af,, the equilibrium interface geometry is like a cylinder. 
Where the axis of the interface coincides with the polar 
axis of the trapping potentials. This occurs when Aj > ctj 
i.e. in the oblate condensates. Here, the phase separation 
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is along radial direction and analogous to planar interface 
in cigar shaped condensates. If p is the radius of the 
interface cylinder, then in TFA 
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Above set of equations define the stationary state of 
TBEC in the oblate shaped condensates. Like in the 
planar geometry, p is the parameter of variation. To ver- 
ify the scheme, we consider pan cake shaped (A^ 3> on 
) TBEC of 85 Rb- 87 Rb mixture in coincident traps with 
Xi as 50.0 and parameter set a. Then from our scheme 
the equilibrium state has cylindrical interface of radius 
5.84a osc . The value from numerical solution of GP equa- 
tion is 5.80a osc . The two results are in very good agree- 
ment and validate our minimization scheme. The con- 
tour plots showing the absolute value of wave functions of 
85 Rb and 87 Rb, obtained by numerically solving Eq.d2j), 
are shown in third image from left in lower panel of FigTTl 



III. ROLE OF THE INTERFACE ENERGY 

In the TFA calculations discussed so far, as mentioned 
earlier, the interface energy is neglected. Accordingly, 
the variational schemes we have adopted incorporate ap- 
propriate interface geometries. However, a general mini- 
mization by considering all the possible interface geome- 
tries favors the ellipsoidal interface as the equilibrium 
configuration. For example, though the cylindrical in- 
terface for A; = 50.0 reproduces the numerical results for 
TBEC of 85 Rb- 87 Rb with parameter set a, the minimiza- 
tion with ellipsoidal interface has lower E. This is evident 
from the values of E, calculated over a wide range of A 
for the three interface geometries, shown in Fig. [3] 

As discussed in ref.|12j. the planer interface observed 
in experiments [3] emerges as the equilibrium geometry 
when the interface energy is considered. When a osc is 
much larger than the interface thickness, the total excess 
energy arising from the finite interface tension |12j is 



n A = 



V(r)] 



3/2 



(14) 




FIG. 3: The upper and lower set of three plots show the 
variation of E with and without the interface energy in TBEC 
of 85 Rb- 87 Rb with parameter set a. Where E of planar (red 
curve), ellipsoidal (black curve) and cylindrical (blue curve) 
geometries are examined as function of A. The inset plots 
show the region around the cross over points. 



Here & are the coherence lengths and K, 
K) are defined as 



K = 
6 = 
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Here the integration is over the interface surface area. 
The above expression is valid provided K > 1.5 and 
C2/C1 < L In the present work we consider TBECs in 
strongly segregated regime with £2/^1 < 1 and hence the 
interface energy in Eq.(14l is applicable. Then the in- 
terpenetration depth is proportional to y/^i/K 1 / 4 and 
— > in the limit 1 /K — » 0. In this limit there is no over- 
lap and TFA solution is an excellent approximation. The 
equilibrium geometry is then the one which minimizes the 
total energy: sum of TFA energy and SI a- 

A precise determination of is essential to ob- 
tain correct geometry of the phase separated TBEC. 
To a very good approximation, the interface en- 
ergy is proportional to interface area. The inter- 
face area in planar and cylin drical geom etries are 
2n((2pi — A 2 L 2 )/a 2 ) and Anp^2pi — a 2 p 2 /Ai respec- 
tively. For prolate and oblate geometries the interface 
areas are 2irR\ + 2TtR\{a\ sin -1 e)/(eAi) and and 2ttR 2 + 
n(a 1 Ri/Xi) 2 ln((l + e)/(l — e))/e respectively. Here the 
cllipticity e are <J\ — (X\ja\) 2 and \/T— («i/Ai) 2 for 
prolate and oblate respectively. The interface areas in 
the three geometries for TBEC of 85 Rb- 87 Rb with pa- 
rameter set a, using our semi analytic scheme developed 
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in previous section, are shown in Fig. [4] The compara- 
tive study reveals that for A <C 1, planar and ellipsoidal 
geometries have lower interface area than the cylindrical 
one. Whereas for A ^> 1, the cylindrical and ellipsoidal 
geometries have lower interface areas. In these two do- 
mains the interface area of one geometry is much lower 
than the other two and hence interface area can decide 
the preferred ground state geometry. For A close to one, 
the difference in the interface areas of the three geome- 
tries is small and surface tension is more crucial than 
interface area to determine the ground state geometry. 
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FIG. 4: Plots showing the interface areas as the functions of 
A, in the TBEC of 85 Rb- 87 Rb mixture, for three geometries: 
planar (red curve), ellipsoidal (black curve) and cylindrical 
(blue curve). The inset plot shows the interface energy. 

In the following subsections we examine the impact 
of £Ia iu two domains: prolate shaped potentials (A^ < 
1 ) and oblate shaped potentials (Aj > 1). For higher 
symmetry and simplified analysis we choose Ai = A 2 = A. 



A. Prolate trapping potentials 

In the A < 1 domain, at low values of A, the ellip- 
soidal geometry has higher ground state energy than the 
planar geometry. As A is increased, keeping the other 
parameters fixed, the ground state energies of both the 
geometries increase. However, the planar geometry has 
higher rate of increase. Then at A Q , which is close to one, 
the energies of the two geometries are equal. Beyond this 
critical value, the energy of ellipsoidal geometry is lower 
and is the ground state geometry. 

For the 85 Rb- 87 Rb mixture with the parameter set a, 
the total energy E and interface energy Qa as functions of 
A for the two geometries are respectively shown in Fig. [3] 
and Fig. [4] (inset plot). Since the value of A a depends on 
the parameters of the system, we examine the variation in 
A a as function of the ratio N%/Ni. For this we fix N\ and 
vary A^. Then calculate A a as a function N2/N1. When 
N2 is decreased A a increases initially and then decreases. 
This is shown in Fig. [5] 
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FIG. 5: Ground state geometry as a function of A c (A a and 
Ab for planar-ellipsoidal and ellipsoidal-cylindrical transition 
respectively) and ratio of population N2/N1 in the TBEC of 
85 Rb- 87 Rb mixture. Inset plots show the variation in energy 
of two lowest energy structures with A for oblate trapping po- 
tentials with interface energy correction. Red and blue curves 
correspond to ellipsoidal and cylindrical geometries respec- 
tively. Each pair curves corresponds to different N2 but same 
Nl The uppermost pair is for Ai = N 2 = 50, 000. The next 
lower pair of curves has A2 = 46, 000, 42, 000 and so on. 



B. Oblate trapping potentials 

In the A > 1 regime, the ellipsoidal or cylindrical in- 
terface geometry is the preferred ground state geometry. 
The planar interface has higher and not favored. For 
A close to 1, the ellipsoidal geometry has lower energy, 
but looses stability as A is increased. This is due to the 
higher rate of increase in the Qa for ellipsoidal geometry 
At the critical value At, and beyond, cylindrical geometry 
has lower total energy and takes over as the ground state 
geometry. 

For the 85 Rb- 87 Rb mixture with the parameter set a, 
the total energy E and interface energy Q,a as functions 
of A for the two geometries are shown respectively in 
Fig. [3] and Fig. [1] (inset plot). For the same parameter 
set a; the value of A&, first decreases and then increases 
on decreasing jVjj. This is shown in Fig. [5] 



IV. CONCLUSIONS 

There are three distinct interface geometries of the 
ground state of TBEC in phase separated regime. We 
have developed a semi-analytic scheme to determine the 
stationary state parameters for each of these and demon- 
strate the validity by comparing with the numerical re- 
sults. We find in TEA, when the interface energy is ne- 
glected, the ellipsoidal geometry has the lowest energy 
for all values of A. Hence is the preferred ground state 
structure. In this structure one species envelopes the 
other and interface geometry and overall shape of the 
TBEC is ellipsoidal. To explain the experimentally re- 
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alized stationary state structures of TBECs we include 
the interface surface tension. We find that minimizing 
total energy, sum of TFA energy and fi^, gives the right 
interface geometry. Then, in our semi-analytic scheme 
the ellipsoidal geometry no longer has the lowest energy 
for all values of A. For cigar shaped traps (A <C 1 ), the 
structure with the planar interface is the ground state 
geometry. While for pan cake shaped traps (A>1) the 
cylindrical interface is the ground state geometry. For 
the values of A close to unity ellipsoidal structure is the 
ground state geometry. 



V. APPENDIX 

In case of planar interface between the two species 
trapped in potentials with separated minima, the expres- 
sions for Ni, iV 2 and E are: 
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These equations reduce to those for coincident centers on 
substituting l\ = L\ and zq = 0. 
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